On Statistics and 1/f Noise of Brownian Motion 
in Boltzmann-Grad Gas and Finite Gas on Torus. II. Finite Gas 
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An attempt is made to compare statistical properties of self-diffusion of particles constituting 
gases in infinite volume and on torus. In this second part, derivation, from BBGKY equations, of 
roughened model of self-diffusion is revised as applied to finite JV-particle gas under micro-canonical 
ensemble. The model confirms existence of characteristic time ~ N, in units of free flight time, for 
cross-over between non-Gaussian and Gaussian regimes of diffusion, but then loses its legacy. 
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I. INTRODUCTION 

In the first part [lj of present work a most sim- 
plified version of the "collisional approximation" Q of 
the Bogolyubov-Born-Green-Kirkwood-Yvon (BBGKY) 
equations was formulated and then used to analyze statis- 
tics of Brownian motion of gas particles in an infinite- 
volume gas under the Boltzmann-Grad limit ( v — > oo , 



0, A 



v-i 



■ const , where T> is space di- 



mension, v and 8 are concentration of gas particles and 
radius of their repulsive interaction, respectively, and A 
their mean free path). Now, our aim is to do something 
similar in respect to finite- volume gas on torus with finite 
total number of particles N . 

Recall that the "collisional approximation" was intro- 
duced in (one can see also [|| or Appendix B in [l[ ) be- 
ing mentioned as correct transition from exact BBGKY 
theory to Boltzmann-like theory in the case of spatially 
nonuniform gas, and in such sense it pretends to real- 
ization of earlier outspoken doubts about validity of the 
Boltzmann equation f4i] or, equivalently, provability of 
Boltzmann's "Stosshalansatz" 5] in this case. One of key 
observations in 0] was that relative motion of colliding 
particles is inner constituent of their collision and there- 
fore should be excluded from its outer characterization. 
In other words, probability density (ensemble- average 
concentration) of pair collisions F 2 undergoes an equa- 
tion dF 2 /dt = -V-dF 2 /dR+... , where V = {v 1 +v 2 )/2 
and R — (ri+r2)/2 are velocity and position of the col- 
lision (more precisely, its center of gravity) , and the dots 
replace "collision operators" acting onto the velocities. 
It is easy to prove that in actually nonuniform situation 
those F 2 (in contrary to general-position pair distribu- 
tion function for non-colliding particles) can not be fac- 
tored into product Fi(t,n ~ R,Vi)Fi(t,r2 ~ R,v 2 ) of 
two one-particle distribution functions. Hence, it turns 
out to be independent supplementary characteristics of 
non-uniform gas. 

If that is the case, then concentration F3 of three- 
particle clusters, or "encounters" @], which enters the 



dots above, also is independent. It in its turn involves 
four-particle encounters, and so on. Thus we have to con- 
sider many coupled linear kinetic equations, which could 
be reduced to the single nonlinear Boltzmann equation 
under strict spatial uniformity only. 

As in (or [3] or [if) we will deal with not "thermody- 
namical" but "statistical" ( "informational" ) spatial non- 
uniformity implied by information about that (at some 
initial time moment t = ) one of gas particles is def- 
initely positioned near some given place r = r . Cor- 
responding (normalized to unit) distribution function of 
the whole gas can be chosen as 



n 



N 



F N (r,v) =F*\r,v)-'£5( v \r j -ro) , (1) 



C N (E) 



N 



9. 
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S(E N (r,v) -E)6W [X>; 1 • 
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where ft — l v is volume of gas contained in "flat cubic" 
torus 0<r Q <^ (a = 1 -j- T>), Em{t, v) is total energy 
of gas, and F^ q (r,v) its equilibrium distribution with 
Cn(E) being normalizing multiplier. Thus, unlike [l|, we 
will rest upon micro-canonical ensemble, with fixed total 
energy and zero total momentum, which looks natural at 
finite number of degrees of freedom. 

We want to study consequent evolution of F^ in or- 
der to obtain statistical characteristics of Brownian mo- 
tion of initially localized particle. Unfortunately, direct 
solution of exact Liouville equation for Fn is impossi- 
ble. From the other hand, the Boltzmann equation (or, 
rather, Boltzmann-Lorentz equation [5|) for F\ can be 
easy solved, but in view of above reasonings (see for 
detail discussion) it seems inadequate with respect to spa- 
tial dependence of F± . Therefore by perforce we have to 
use the theory from |l[ adapting it to finite system. At 
that, we can not exploit the Boltzmann-Grad limit which 
helps to obviate most of difficulties of the "collisional 
approximation" descending from vagueness of the very 
concept of "collision" Q. By these reasons, at present 
only quite imperfect and caricature collisional model of 
nonuniform finite gas is available. Nevertheless, may be, 
it is better than uncritical use of the Boltzmann equation, 
all the more, better than nothing. 
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II. MODEL EQUATIONS 

1. Let us make non-principal assumptions as follow. 
First, of course, that total number of particles is large 
enough, N 3> 1 . Second, mean free path is comparable 
with the torus dimensions: A ~ l/uS 13-1 ~ I , with v = 
N/fl = N/l v . Under these assumptions, estimate of a 
number of simultaneously occurring collisions obviously 
yields a value much smaller than unit: N(S D N/Vt)/2 ~ 
jSf-i/CD-i) ^ ^ _ Hence, our gas is sufficiently rarefied. 
On this ground, when considering probabilities of particle 
velocities, we can neglect potential energy contribution 
to E N (r, v) and write in E = NVMv%/2 , with v a 
being characteristic velocity (M is particle mass) . 

In order to get ability of watching for full rotations 
of particles around torus, let us use extended configu- 
rational space — oo < r a < oo , so that its cubic cell 
(m a — 1/2)1 < r a < (m a + 1/2)1, with integer numbers 
m a , corresponds to m a rotations along a-th dimension 
during time period after t — . Under this represen- 
tation of life on torus, potential energy of interaction of 
any two particles changes to periodic function of coordi- 
nates, U(p(ri — rj)), where p(r) means r modulo I: 
Pa(r) = r a —m a l when (m a — 1/2)1 < r a < (m a + l/2)l . 
At that, distribution ([T]) also will be extended to the 
whole space — oo < r a < oo , with Q now being men- 
tioned as its (infinite) volume. Such statistical ensemble 
implies, clearly, that rotations of only a sole (initially 
localized) particle will be watched, while getting of in- 
formation about absolute positions and displacements of 
other particles is beforehand eliminated. 

Next, consider particular n-particle distribution func- 
tions which follow from |T]): 

F n (n...v„) = ... F N (r,v) = 

Jn + l JN 



~ N 

(N - n) 
N 



F^( ri ...v n )J2^ V Hr,-r Q ) + (2) 

.7=1 



F n l 1 (ri...r n ,r ,vi...v n ,v )dv , 



) , (3) 



F^(n...v n ) = [ ... [ F£(t,v 

Jn + l JN 



where f ... — J...dv s dr s . Since, according to fT}, 
F n ~ , and il has become infinite, while is finite, 
it is evident that the second term on right-hand side of 
([2]) factually vanishes in comparison with first term and 
should be thrown away when considering nonuniform (lo- 
calized) parts of FnS. Besides, it is reasonable to move 
off also common multiplier 1/N . Thus we can write 

n 

F n ( ri ...v n ) = ilF?(r l ...v n ) t {V) (r 3 -r Q ) (4) 
with normalization J ...J F n — n£l n ~ l . 



2. Formula ((4]), along with |T]) and ([3]), presents initial 
conditions for the BBGKY equations 



dF n 
dt 



/n 
Y^Ljn+lFn+l , (5) 
-+1 i= i 



9r, 



1< i<j <n 

_d d_ 

dp j dpi 



In "collisional approximation" we want to transform the 
right-hand integrals in ([5]) into Boltzmann-like "colli- 
sional integrals" . 

This is possible in those cases only when statistical 
ensemble under consideration contains all various con- 
secutive phases (instant dynamic states) of a collision 
process, moreover, when all that phases are represented 
by the ensemble with equal probability. Otherwise, the 
collision process can not be replaced by momentary col- 
lision event, since the latter must conserve probability 
(in essence, conserve particles). Of course, hardly a 
non-trivial ensemble ever (even at Bogolyubov's "kinetic 
stage" ) undergoes such requirement in formally rigorous 
sense. Therefore this requirement should be interpreted 
as instruction how to make "collisional approximation" . 

Firstly, we must consider special values of F n ( n > 1 ) 
corresponding to clusters ("encounters" 0, 0]) of n mu- 
tually close particles (e.g. satisfying, in real torus space, 
conditions \n — Tj | < d , where d is significantly smaller 
than / but greater than 6 ) . At n = 2 the matter con- 
cerns either literally collision or close encounter of two 
particles with no essential interaction, while at n > 2 
we speak about chains of close pair collisions or their 
combinations with mere encounters (in particular, might- 
have-been or "virtual" collisions in the sense mentioned 
in i or g). 

Secondly, divide term L^ n 'F n in ([5]) into two parts, 
one coming from non-uniformity of statistical ensemble, 



that is FnS dependence on translations Tj 



Ar , 



and other from its dependence on dynamic phase of n- 
particle encounter. This can be done in the form 



L^F n = -vw ■ y ^ - 

. or 4 



=i 



dFn 



(6) 



where V^ n ' = (t>i + ... + v n )/n is velocity of centroid 
of the n-particle cluster, and O is its inner time in the 
centroid's reference frame. Obviously, ^d/drj is just 
infinitesimal operator of plane-parallel displacement of 
the whole cluster. About O see Q or Appendix B in 

Thirdly, and chiefly, following the above requirement, 
second term in © must pass for zero: dF n /dQ — > 0, 
since the concept of collisions imposes equiprobability of 
its dynamic phases (in particular, its border "in-" and 
"out-" states), i.e. conservation of probability. 

After all, F n , being taken at |/o(t", — Tj)| < d, be- 
comes probabilistic characteristics of n-particle encoun- 
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ters as one-piece events. Two right-hand terms in §5§ 
turn into drift and source of the encounters, respectively 
(to be precise, drift and source of their ensemble-average 
concentration). At that, the requirement dF n /8Q — > 
once again plays important role: it serves for transform- 
ing the source into sum of collision operators which act 
upon F n+1 

However, still F n possesses too many arguments (after 
removal of , only by unit less than generally) , including 
impact parameters and other internal geometric details 
of an encounter. Therefore simultaneously let us exclude 
the latter arguments by averaging F n over them. Before- 
hand notice that spatial dependence of exact solution to 
([5|) starting from (j4|) must look like 

F n = F n (t,n; p(r 2 -r 1 ),...p(r n -r 1 ))+ CP , (7) 

where F n is symmetric with respect to indices 2 n and 
"CP" means sum of n — 1 cyclic permutations of all n 
indices. For this reason, it is sufficient (with no loss of in- 
formation) to average F n over only r.j belonging to one 
and the same cell of the extended space (and multiply re- 
sults by some constants to keep suitable normalization). 

Thus, fourthly, introduce W(r) as some spherically 
symmetric bell-shaped (and normalized to unit) weight- 
ing function with width ~ d (i.e. approximately enclos- 
ing a region assigned to collisions and encounters) and 
then averaged distributions as 



essence this is not significant defect of our model, be- 
cause, as we will see, factual role of collisions is exchange 
between momentum of "outer (n + l)-th" particle and 
total momentum of n-particle cluster. 

What is more important, we came to the point where, 
fifthly, we need in a hypothesis like Boltzmann's Stossha- 
lansatz ("molecular chaos") in order to reduce func- 
tions F m to functions F and thus close equations 
(|9"|). But both trivial reasoning mentioned in the In- 
troduction and deeper considerations 0, 0] show that it 
can not be literally Boltzmann's hypothesis. If speak 
about infinite gas, nothing prevents a cluster's ("j-th") 
particle and exterior ("(n + l)-th") particle be sta- 
tistically independent in the sense of their velocities, 
but, naturally, their positions are statistically depen- 
dent, which is expressed by F n+1 (t, R, V\ ... v n , i>„+i) = 
F (v n+1 )J F n+1 (t,R,v 1 ...v n ,v')dv' (with F (v) being 
Maxwell distribution when gas is in thermodynamical 
equilibrium). However, in finite gas, framed by the 
micro-canonical ensemble (JXJ) , even velocities of "j-th" 
and "(n + l)-th" particle must be statistically dependent. 

In view of this new difficulty, relations of F^a-i to F n+ i 
will be analyzed in concrete way of further simplification 
of the model. 

3. According to JT]), (fJJ and ([H]) and notions made in 
the beginning of this Section, initial conditions to equa- 
tions (T9l) must be 



(l n - 



F 



l[W( 

j'=i 



Tj — R) drj 



(8) 



Applying such operation to (O, in view of (O, we can 
replace expressions J dr n+1 U{p(rj — r„ +i )) ... by ex- 
pressions (v£l/N) J dr n+ i U(r 3 ■— r n +i) ... with f2 being 
(infinite) volume of the extended space and v (real finite) 
concentration of gas. Taking into account also equalities 
([6]) and dF n /d& = , and manipulating by analogy with 
it is not hard to obtain 



d_F\ 
dt 



dF n 



N - 



dR 



N 



Sjn+\F n+1 (9) 



with y(") = n~ l Y?i 



Symbol Sj n +i designates 



Boltzmannian linear collision operator describing colli- 
sion between j-th member of n-particle cluster and an 
exterior (n+ l)-th particle. The latter came from outside 
of the cluster and is in incoming state with respect to it, 

being placed somewhere at its border. Function F^l_ 1 
just represents such situation. 

The last collision terms in ([9]) split to purely pair colli- 
sions (without participation of other members of a clus- 
ter at n > 2 ). Strictly speaking, this statement is simply 
demonstrable only for infinite gas under the Boltzmann- 
Grad limit (when d can be chosen so that S/d — * 
despite d/X — > and any cluster becomes "infinitely rar- 
efied from within" although "infinitely dense from with- 
out"). Nevertheless, even if it is formally wrong, in 



F n (t = 0,R,v 1 ...v n ) = nW{R)F e n q {v 1 ...v n ) , (10) 
F e n q (vi...v n ) = J dv n+1 ... J dv N x 



J'=l 



vi=i 



where E = NT>Mv%/2, we took r = , and F^ 
are equilibrium velocity distributions normalized to unit. 
The integration yields 



F„ oc 



(gf=i V; 
N-n 



(N-n-l)T>/2- 1 



Soon we will use the first and second statistical moments 
of these distributions: 



(Vj) = , (viaVjp) = Vq 5, 



a/3 



N — 1 



, (11) 



where ([...]) = / — /[—] F^ dv\... dv n and i,j < n, 
and besides the incomplete integral 



v n +i F^ q +1 dv n+ i 



Vl 



N-n 



V n —eg 



(12) 



Obviously, factor — (vi + ... +v n )/(N— n) here is nothing 
but conditional average value of i> n +i under fixed Vj< n . 
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4. As in [l|, mostly we are interested in spatial distri- 
butions and conjugated probability flows 



W n (t, R) 



F n (t, R, v\...v n ) dvi... dv n 



J n (t,R) = J ...J F n {t,R,v x ...v n )dv 1 ...dv n , 

first of all, Wi(t, R) which is probability density of ran- 
dom Brownian displacement of initially localized particle 
(from its startposition near R = ) . 

Similar to let us derive from © roughened but 
closed and solvable equations for these distributions, by 
using the fact that after a few collisions (counted off 
from t = ) cross-correlation between previously accu- 
mulated displacement R of Brownian particle and its 
current velocity v becomes small. Indeed, the condi- 
tional average value of v a at given R is ~ R a /t . Since 
VQTft (with D being diffusivity and 
Tf ~ X/vo mean free flight time), modulus of this aver- 
age is ~ vq y/2rf/t , hence already at t ~ 20t/ it is small 
enough in comparison with magnitude vq of equilibrium 
velocity fluctuations. 

At larger time the displacement-velocity cross correla- 
tion can be treated as a weak perturbation. Clearly, cor- 
responding weakly perturbed distribution functions F n 
can be written as 

F n (t,R,Vi...v n ) = Fl q {v x ...v n ) x (13) 
x [W n {t,R) + w n (t,R) ■(«! + ... +v n ) + ...} , 

where first term in squire brackets represents zero-order 
contribution from asymptotic with independent velocity 
and displacement, while the dots replace higher-order 
contributions. The latter will be neglected and omitted. 
Substituting (fT3|) into definition of the flows and using 
(fTTj) we find relations of w n to them: 



J n [t, R) 



N ■ 



N- 1 



Vq w n (t, R) 



(14) 



Next notice that the question about connections be- 
tween distribution functions F^ 1 and F n+ i will be 
simply resolved in the framework of the first-order ap- 
proximation if we write 

F^. 1 (t,R,v 1 ...v n+1 ) = F^ l q +1 (v 1 ...v n+ i) x (15) 

n 

x { W n+1 {t, R) + w n+1 (t, R) ■ a n Vj + c n v n+1 }} 

3=1 

and choose coefficient a n and c„ to be such that mean 
value of D n+ i turns into zero, 



v n+1 F n+1 dv 1 ...dv n+1 = 



(16) 



while mean values of Vj with j < n + 1 remain exactly 
the same as under distribution F n+ \ : 



vj F n+1 dvi... dv n+1 = / vj F n+ i dvi... dv n+1 



With the help of (jTT|) - (fT2|) . these requirements yield 
n N -1 



N — 1 ' 



N 



(17) 



respectively. 

At such a n and c„ , expressions (fT5l) and p6|) say 
that velocity of outer (n + l)-th particle is not correlated 
with position of n-particle cluster it came in (as if it 
was a particle from equilibrium thermostat), although 
its position can be correlated (as if it was full member 
of greater (n + l)-particle cluster). This is just what is 
required by the "weakened molecular chaos" suggested in 
0, HI as a compromise of mutually contradicting "pure 
Stosshalansatz" and conservation of probabilities. 

Consider time derivatives of the probability flows, com- 
bining |9]) with expressions (113|) and (fl~5| . noticing that 



N-n 
n(N-l) 



(18) 



due to (|lip . and everywhere keeping only lowest-order 
nonzero contributions: 



dJ n 



N-n 



dW n 



dt n{N- 1) dR 



N 



dvi... / dv n V {n) 



(19) 



x w n+1 (t,R) ■ Si n +i [ ^2 Vj + c n v n+1 ] F* q +1 
Action of the collision operator is defined as usually, 



'in+l 



db / dv n+ i x 



X \Vj — V 



n+1 



[/(... v[. 



J n+1 



)-/(• 



,Vi...V 



n+1 



)] 



where b is (T> — l)-dimensional impact parameter, and 
v[ = v' l (v l ,v n+1 ,b) and v' n+1 = v' n+1 (v l ,v n+1 ,b) are 
those in-state velocities which result in given out-state ve- 
locities Vi and v n+ i after collision. Taking into account 
that Vi + v n+ i = v[ + v' n+1 and F n+1 are invariants of 
collision, it is easy to verify that 



3=1 



S in+l 

db I dv n+ i[ 1 — cos 9}\vi 



1 



Vj + c n v n+1 ] F n+1 =-~(l- Cn )vx 



J n+1\ 



\~= eq 
Vn+l)-f n+1 



where 9 = 0(\Vi — v n+ i\ ,b) denotes scattering angle as 
a function of collision parameters. Of course, spherical 
symmetry of interaction potential was assumed. Now 
perform in (fT9|) integrations first over all Vj with j ^ 
i,n+ 1 , using relations similar to (|12|) . and second over 
Vi and v n +i , paying attention to symmetry of resulting 
equations, and introduce the "relaxation rate" 



(N- l)v 



ANvl 



dv / dv' 



db \v-v'\ 3 



[l-cos9(\v-v'\,b)]F° q (v,v') 
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Finally, with the help of (14]) and (171), equation l[l9l 
transforms to 



dJ„ 



v 2 N -n dW n N -n 



7- 



9* n JV-1 9iJ ' N 

It remains to supplement (|20p by equations 

9* 



, (20) 



(21) 



OR ' 
and initial conditions 

W n {t = Q,R) = nW(R) , J n (t = 0,H)=0 , (22) 



which directly follow from (9|) and (fTO)) . 

Clearly, 7 can be identified with inverse mean free 
flight time 1/t/. We may replace VF(i2) by delta- 
function 8 V (R) , since anyway we have to consider spa- 
tial scales much greater then d . Of course, index n in 
(l2"0")) - (|22p takes values 1 < n < N — 1 , because zero total 
momentum in our micro-canonical statistical ensemble 
(T|), (p~0|) requires to set 



J N (t,R) = 



(23) 



The same is dictated by (H]), or (H at n = N . 

If we chose canonical ensemble, then in (|16l) we would 
have to write c„ = and a n = 1 in place of p7[) . while 
formulas (|14[) and (T5|) would change to J„ = UqIO„ and 
(Vq ( ™V /3 ( ™ ) ) = <5 Q/3 Uo/ n i resulting in 



dt 



v 2 dW n N-n 



7 . 



J n 



+ 1 



(24) 



n OR ' N 
instead of PU|) . now with 1 < n < N and «7/v 7^ . 

III. MODEL ANALYSIS 

1. First, as in 1], consider solution to equations 
(|2"2")l in terms of Fourier and Laplace transforms 

E n (t, k) = J cxp(ifc ■ r)W n (t,r)dr , 

$„(p, k) = / exp(-pt) 3 n (i, fe)dt 
Jo 



At that, in view of obvious isotropy of Brownian motion 
under interest, it is sufficient to watch for its projection 
onto some axis only, formally changing "wave vector" k 
by scalar k . Then 

1 00 (—k 2 \ n r°° 

with (R 2n (t)) = J R 2n W 1 (t,R)dR being statistical mo- 
ments of (a projection of) Brownian displacement. 



Now from (I20l)-(I22I) we come to finite series: 



/ ,x v 2 k 2 N , . 

p$ n (p, k)-n= x 25 

pj(N- 1) 



AT-1 s 

><e n 



(AT - m)7 



Np 



v 2 k 2 (N - m) 
mp 2 {N - 1) 



Correspondingly, Laplace transform of (R 2n (t)) is (2n + 
./V— l)-order polynomial of 1/p with powers from 2n + l 
to 2n+JV-l. Hence, (R 2n (t)) itself is (2n+7V-2)-order 
polynomial of t with powers from 2n to 2n + A^ — 2 . 
In particular, 



S=l 



_jlV 
at p ; 



which means that 
d 



dt 



(R 2 (t)) = 2D 



1 



1 



21 

N 



N-l 



(26) 



with diffusivity D presented by 
D = 



v 2 N 



7 (iV-l) 



Consider also fourth-order moment of (projection of) 
the displacement: 



e- pt (i? 4 (t))^ 



24L> 2 7 



p" 



N-l 



V (N-l)\ f 7_V 
£t (iV-s-1)! I AW 

which yields in the time domain 



s 1 

E 1 

^ — ' rr. 



m— 1 



S 

m AT 



-^<i? 4 (t))= 24Z) 2 ( tf^ 1 



1 N — 1 



N 



A^ 



AT-1 



24D 2 ^ 



JV 



N-s 



(27) 



with factor d = defined by 
1? = 1 



7t 
A^ 



Polynomials representing higher-order statistical mo- 
ments may be considered in terms of functions Vn m \t) 
introduced by analogy with right-hand sides of (f2"6"l) - ([2"T)) : 



2m 



dt 



m — 1 



J R 2m ^J n (t,R) dR 



(28) 



dt 



R 2m W n (t,R)dR ee (2m)! D m V^ m \t) 



6 



CM 

Q 

H — 1 

DC 
T3 




100 



FIG. 1: Plots of (24D 2 )~ 1 d 2 (R i (t))/dt 2 via -ft correspond- 
ing to l|27p at N = 21 (blue) and N = 51 (black) in compari- 
son with In (-ft) + C (red) which approximates limit N — ► oo 
( C is Euler constant). 



From ([20 |) -(|22 p recurrent equations for Vn (t) directly 
follow: 



dvir ] 

dt 



N -n 



7- 



N 



I T/(m-l) _ t/M 



(29) 



with 1 < n < JV - 1 , K (0) (t) = n , m > 1 , initial 
conditions V^ m ' > (t = 0) = , and Vy" 1 ' = because 
of (|2"3")) . In particular, at m = 1 equations (|2"5|) yield 
generalization of (|26[) . 



^ /" R 2 dR = 2D [1 







N-nl 



(30) 



2. Expressions (I26p and (|27p look quite formally satis- 
factory and physically reasonable as long as -ft < N , that 
is > . Moreover, as it naturally might be expected, 
at sufficiently large N and 1 <C -ft < N (more precisely, 
at \$\ N < 1) they approximately coincide with expres- 
sions for infinite system [l[ , that is 2D and 24D 2 In -ft , 
respectively (as the consequence, while -ft < N the dif- 
fusivity l//-noise exists in the sense explained in [H, HJ). 
Similar statements are true also in respect to various 
higher moments J R 2m W n (t, R) dR which follow from 
([23)) or (f2"9"| and to the whole probability distribution 
functions W n (t, R) (at least at N — n ^> 1 ). Particular 
example is shown in Fig.l. 

However, at -ft > N ( d < ) the same expressions 
(|2"tj)) . ([77)1 and (|30p become obviously senseless. Notice 
that $ and its various powers $ s play roles of relaxation 
factors being direct analogues, especially at s ~ N , of 
the exponent exp (— -ft) in infinite system. Therefore, 
negative $ with growing absolute value has none physical 
meaning. 

Hence, we should conclude that solution to equations of 
our model (|20p - (p22")) possesses peculiar time point N/-f = 
Nrf . Out of it the model does not work. 



These fact can not be addressed to some incidental 
details of the model. Indeed, other model correspond- 
ing to canonical ensemble and represented by equations 
in place of (|2"01) possesses the same peculiar point. 
It is not hard to trace that the peculiarity arises from 
two things: presence of factor N — n in collision term of 
equations (|20p and characteristic open-chain structure of 
these equations with n-th member referring to (n+l)-th. 
Both have migrated from the basic BBGKY equations. 
But the second is also child of the "weakened molecular 
chaos" hypothesis. Hence, just the latter (to be precise, 
its concrete realization on the form of (fT5|) - ([!())) ) is re- 
sponsible for the peculiarity. 

It seems natural to interpret this situation as evidence 
that at -ft = N period of non-Gaussian behavior of our 
random walk is over. Then predicted peculiar time N/-f 
conceptually as well as numerically coincides with time 
T ng suggested in |l[ as characteristic time scale separat- 
ing "non-ergodic" (non-Gaussian) and "ergodic" (asymp- 
totically Gaussian) stages of the walk. 

3. In order to construct a better model (where weak- 
ened molecular chaos would self-destruct with time into 
complete one) we should return to higher level of equa- 
tions ([U or BBGKY equations. This is subject of sepa- 
rate investigation (principally, it is plausible if "trivial" 
Gaussian asymptotic creates nontrivial problems). 

At present, we confine ourselves by discussion of most 
weak spot of the model. Undoubtedly, that is the "end 
cap" (|23p which makes closure of the chain (|2l7)) . but 
rather bad closure, because it means that (N—l) -particle 
cluster loses any interaction with rest iV-th particle. 
Here is the reason to remove (N— l)-th of equations (|20|) 
from affected zone of the "weakened molecular chaos" hy- 
pothesis (|16p - (fT7|) . Then it takes the form 



dJ 



N-l 



r 2 
' 



dW N -i 



1 



01 



{N-l) 2 dR 
= v 2 a N -i(l - c N -i)w N 



(31) 



Equality (|23p remains true, but now it is out of work. 
The model becomes temporarily non-closed, and we get 
chance to correct its behavior after t ~ r ng (notice that 

r ng is just time necessary to transfer influence by Jn 
from last to first of equations (|2"0"|) ). 

Let us demonstrate that in presence of Jn the model 
comprises Gaussian asymptotic at t/T ng — > oo , more- 
over, without significant assumptions about accompa- 
nying asymptotic of Jn ■ According to ([2"5]). Gaussian 
asymptotic means that v} m \t — > oo) — > 1 . Then other 
variables in equations (f2"9"|) also have fixed points fully 
determined by (oo) (as combined with V^ ' = n): 

( ) f (n-m)!/(n-l)! , 1 < m < n 

V n (°°) = i T ,(m-n+l)/ w, n , 

L Vi (oo)/(n — 1)! , m> n 



(32) 

Here 1 < n < N , and quantities (t) repre- 

sent J]y(t,R) by the common rule (|28p with n = N 
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and Jpf(t,R) in place of Jpf(t,R). We did not ful- 
fil substitution V^ m ^(oo) = 1. This helps to visualize 
that \ r k n \oo) at least at m < N appear insensible to 

V-[ (oo) (i.e. independent on assumption about gaus- 
sianity of the asymptotic) thus justifying above state- 
ment. Notice also that equalities (oo) = 1 imply 

equalities ^"^(oo) = Vi m \oo) (at m > 0), which mean 
that asymptotically J2(t,R) — > J\(t,R) and therefore 
first pairs of equations (f20|) - (|2Tj) stands apart, as it must 
be in standard kinetics. 

IV. CONCLUSION 

In the first part of this paper [l| a simple kinetic model 
of self-diffusion (Brownian motion) in infinite-volume 
gas was developed. It exploits collisional description of 
interaction between particles but declines Boltzmann's 
"molecular chaos" hypothesis because the latter appears 
incompatible with two doubtless theoretical requirements 
to be fulfilled: probability conservation during collisions 
and necessity to deal with spatially nonuniform statisti- 
cal ensemble when considering random walks of gas parti- 
cles. Instead, the "weakened molecular chaos" hypothesis 
suggested in 0] was used, which is consistent with both 
the requirements and claims statistical independency of 
colliding particles in respect to their velocities only but 
not their coordinates. It resulted in a chain of kinetic 
equations whose solution shows that statistics of random 
walks of gas particles does not obey the law of large num- 
bers and remains essentially non-Gaussian at arbitrary 
long time scales, including what can be named 1/f fluc- 
tuations in diffusivity. 

Physical origin of so wild statistical freedom is absence 
of actual cause- and- effect correlations between successive 
fragments of the random walk, in spite of statistical long- 



living correlations formally describing the same freedom 
0, y]. Its mathematical origin is absence of relaxation 
(in the ordinary sense) terms in kinetic equations. Usual 
place of relaxation terms is occupied by references to next 
higher-order equations. This chain of references mean 
that a particle whose random walk is under observation 
constantly accumulates (actual) correlations with more 
and more new particles. 

In infinite gas this process lasts endlessly. In finite TV- 
particle gas it ends at time ~ N (counted from start of 
the observation, in units of mean free flight time) when 
collisions of observed particle with all others have re- 
alized, therefore, initial conditions of the whole system 
have snapped into action. Further motion of this particle 
is fully predestined by already observed walk, represent- 
ing its reflections (mappings) in complete non-decaying 
actual correlation with all the gas. This creates the nec- 
essary prerequisites for decay of statistical correlations 
between particles and realization of the law of large num- 
bers and Gaussian statistics at next longer time scales. 

The model investigated in this second part of the pa- 
per certainly detects time w N of this transition but, 
unfortunately, is unfit for description of the next Gaus- 
sian asymptotic of the random walk. However, in prin- 
ciple, the exact BBGKY equations give a firm base for 
proper improving the model. From the other hand, prac- 
tically sooner just the former non-Gaussian stage is of 
interest, since N means rather large time. In this re- 
spect, the model needs in improvements too. Proba- 
bly, at present form it somehow overestimates "degree 
of non-gaussianity" of the random walk. The early phe- 
nomenological construction of "quasi-Gaussian" walk Q 
(see also references therein and Q) produced much softer 
non-gaussianity (although also with 1 /f noise in diffusiv- 
ity). In any case, there are many interesting questions 
for future. 
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